Triviality of the Ground State Structure in Ising Spin Glasses 
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We investigate the ground state structure of the three-dimensional Ising spin glass in zero field by 
determining how the ground state changes in a fixed finite block far from the boundaries when the 
boundary conditions are changed. We find that the probability of a change in the block ground state 
configuration tends to zero as the system size tends to infinity. This indicates a trivial ground state 
structure, as predicted by the droplet theory. Similar results are also obtained in two dimensions. 
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Controversy remains over the nature of ordering in 
spin glasses below the transition temperature, Tc, and 
two scenarios have been extensively discussed. In the 
"droplet model" proposed by Fisher and Husetl (see also 
Refs. H-gl), the structure of "pure states" is predicted to 
be trivial. This means that there is a unique stateH in the 
sense that correlations of the spins in a region far from 
the boundaries are independent of the boundary condi- 
tions imposed. As a-|C|Qnsequence, the order parameter 
distribution functionOlH, P{q), is also trivial, i.e. is a 
pair of delta functions at q = ^Qea where qea is the 
Edwards-Anderson order parameter. In the alternative 
apprOjach, one assumes that the basic structure of the 
ParisiErLI solution of the infinite range model applies also 
to realistic short range systems. In this picture, P{q) 
is a non-trivial function because many thermodynamic 
states contribute to the partition function, i.e. the pure 
state structure is non-trivial. Monte C^ijl£| simulations 
on short range models on small Iatticeal3il3, find a non- 
trivial P{q) with a weight at q = which is independent 
of system size (for the range of sizes studied) , as predicted 
by the Parisi theory. 

Most numerical work has concentrated on P{q). By 
contrast, here we attempt to determine the pure state 
structure by investigatiag whether spin correlation func- 
tions in a finite regionu far from the boundary, change 
when the boundary conditions are changed. It is inter- 
esting to investigate this question even at T = 0, where 
there are efficient algorithms for determining ground 
states, even though P{q) is trivial in this limit (for a con- 
tinuous bond distribution). Here we show that the Ising 
spin glass in three^ dimensions which has a finite transi- 
tion temperaturetJiLJ Tc, has a trivial ground state struc- 
ture. We also find a trivial ground state structure in the 
two-dimensional Ising spin glass which has a transition at 
zero temperature with long-range order at T = 0. Some 
of our results in two dimensions have also been reported 
elsewherelia (referred to as PY). Similar results for two 
dimensions, as well as results for some three-dimensional 
models (but none for a spin glass with a finite Tc) have 
also been found by Middletonta. 

The Hamiltonian is given by 



where the sites i lie on a simple cubic {d — 3) or square 
lattice {d = 2) with N = L'^ sites (L < 10 in 3d, L < 30 
in 2d) , Si — ±1, and the Jy are nearest-neighbor inter- 
actions chosen according to a Gaussian distribution with 
zero mean and standard deviation unity. We determine 
the energy and spin configuration of the ground state 
for a given set of bonds, initially for periodic boundary 
conditions denoted by "P" . Next we impose anti-periodic 
conditions ("AP") along one direction, which is equiva- 
lent to keeping periodic boundary conditions and chang- 
ing the sign of the interactions along this boundary, and 
recompute the ground state. Finally we change the sign 
of half the bonds at random along this boundary, which 
we denote by "R" . 

To determine the ground state in three dimensioi3S..we 
use a hybrid genetic algorithm introduced by PalElEJ. 
Starting from a population of random configurations 
("parents"), new configurations ("offspring") are gener- 
ated by recombination (triadic crossover) and mutation. 
The population is progressively reduced, with a bias to- 
wards the offspring with lower energy. The algorithm is 
hybridized with a local optimization of the offspring (see 
Ref. EO for details). For each sample and boundary con- 
dition, we repeat the algorithm n^ times, see Table I, and 
take the lowest energy state found. Our values for the 
average ground state energy are in agreement with those 
of Ref. ^. For L < 6, we checked all our results with a 
different method, which consists in repeating many times 
the microcanonical simulated annealing algorithm intro- 
duced in Ref. E3. This algorithm was also used to prepare 
the initial populations of the genetic algorithm for L > 8. 
We discuss later additional checks that we performed for 
the largest sizes L = 8 and 10. 

In two dimensions, we used the Cologne spin glass 
servers, which calculates exact ground states of the Ising 
spin glass with periodic boundary conditions. 

In order to study the dependence of the spin configura- 
tion on boundary conditions we consider a central block 
containing Nb — Lg spins. We compute the block spin 
overlap distribution P^a{q), where a and f3 denote two 
boundary conditions, P, AP or R here, and 
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JE) 
-106.59(4) 
-210.26(7) 
-364.9(1) 
-868.1(2) 
-1697.5(4) 



(E) Ref. ^ 
-106.609(9) 
-210.22(3) 
-364.89(5) 
-868.1(2) 
-1698.8(8) 



TABLE I. Parameters of the simulation and results for the 
ground state energy in d = 3. As a function of the size L we 
show: number of samples Ns, number of runs per sample Ur, 
and the avpage ground state energy (E), from our data and 
that of PaO. 
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is the overlap between the block configurations with a 
and /3 boundary conditions, S" is the value of Si in the 
ground state with the a boundary condition, and the 
brackets (• • •) refer to an average over the disorder. 

Since we work at T = 0, each sample and pair a, (3 
gives a single value for q. The self overlap distribution, 
P^(q), has weight only at q = ±1, since the ground 
state is unique for a given boundary condition. P^{q) is 
normalized to unity i.e. J P^{q) dg = 1, it is symmetric, 
and the allowed g-values are discrete with a separation 
of of Aq = 2/N^, so Pi'„(±l) = iV-S/4. 

If the configuration in the block changes when the 
boundary conditions are changed from a to /3, then 
the block overlap, g^^, will no longer be ±1. Hence 

1 - Papil) / Paai^) IS the probabifity that the block 
ground state changes on changing the boundary condi- 
tions. We will see that 
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with A > 0, showing that the block ground state config- 
uration is uncbaBged in the thermodynamic limit. Our 
interpretationll3lij of this result is that the boundary con- 
dition change induces a domain wall of fractal dimension 
df = d — X, and L~^ is then the probability that the 
domain wall intersects the block. 

Now we discuss our numerical results, for which ake 
Lb = 2. First of all. Fig. |l| shows the root mean square 
energy difference, AE, between P and AP and between 
P and R boundary conditions in d = 3. Apart from the 
smallest size, L = 3, the data for AEp^ap, are consistent 
with the power law variation 
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where the asymmetric error bar comes from systematic 
effects discussed below. The positive value of 9 shows 
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FIG. 1. Upper figure: a plot of the root mean square 
ground state energy differences IS.Ep^ap and I^Ep^r for dif- 
ferent sizes up to 1/ = 10 in d = 3. The fit for the P-AP data 
omits the L = J, point. Lower figure: a scaling plot of the 
symmetrized distribution of the values for SEp^ap, the en- 
ergy difference between periodic and anti-periodic boundary 
conditions for a single sample. The distribution scales quite 
well with the value of 6 obtained from the fit to AEp^ap- 



that the system is stable against breaking up into large 
domains of little energy, which-joiplies that Tc > 0, 
in agreement with earlier workEj Ell. r^jhe value for 9 
is a little larger than earlier estimatesEa for the Gaus- 
sian distribution considered here, but these calculations 
used a much smalleE-range of sizes. For the ±J dis- 
tribution, HartmannE3 studied sizes up to L = 10 and 
found 9 ~ 0.19 ± 0.02, which is just consistent with the 
value here. It is expectedEJ that the distribution of energy 
differences, SEp^ap has the scaling form P{5Ep,ap) — 
L~^ P{6Ep,AP / L^) and this works quite well as shown 
in the lower part of Fig. 1^. As discussed in PY, AEp^p 
is expected to vary as jj-'^~^)l'^ , (= L here), for large 
L. Fig. n^ shows that our results are consistent with this 
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FIG. 2. A plot of histograms of block spin overlaps for 
L = 4 and 10 in d = 3, with block size Lb = 2. Note that 
the allowed values of q are 0, ±0.25, ±0.5, ±0.75 and ±1. The 
left hand column is for the P-AP overlap and the right hand 
column for the P-R overlap. The top row is for L — 4 and 
the bottom row for L = 10. The data is symmetrized and 
normalized so that the area under the histograms is unity. 

Some representative histograms of the block overlap 
distributions are shown in Fig. g for d = 3. One sees that 
for both P-AP and P-R overlaps the weight a.t q — ±1 
increases with increasing L. 

The probability that the configuration of the block 
changes when the boundary conditions are changed, see 
Eq. dj), is plotted for a range of sizes in the upper part 
of Fig. ||. For comparison, the lower part of the figure 
shows similar data in d = 2 (see PY for related results). 

In d = 3, the data for P-AP overlaps lies on a good 
straight line with a slope —A where A = 0.32 ±0.02. This 
leads to a fractal dimension of the domain walls given by 



df =d- A = 2.68 ±0.02 



(7) 



The data for P-R overlaps for larger sizes lies parallel to 
this, as shown in the figure, but with some deviations for 
smaller sizes. However, given the statistical uncertain- 
ties, one can also fit all the P-R data to a power law with 
a different slope giving A = 0.23 ± 0.02. If the ground 
state structure in d = 3 were non-trivial the data for 
the P-R boundary conditions in the upper part of Fig. |^ 
would eventually saturate at a finite value for L — > oo. 
There is no sign of such a saturation for the sizes that 
we are able to study. We also considered fits of the form 
a + bL^^ , finding that, although a range of positive val- 
ues of a is not ruled out by the data, the value a = is 
statistically preferred. 
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FIG. 3. The figures show the probability that the spin 
configuration of the central block changes when the bound- 
ary conditions are changed. Upper figure: Data for three 
dimensions, for different sizes up to L = 10 with block size 
Lb equal to 2. The solid line through the P-R data is parallel 
to the fit for the P-AP points (slope = —0.32), and fits the 
results for large L. The dashed line is fit to all the P-R data 
and has a different slope of —0.23. Lower figure: a similar 
plot but for d = 2, with sizes up to L = 30 and L^ — 2. 



We also investigated the effects of changing the bonds 
on one boundary to other random values (rather than 
just changing the sign) and also changed the bonds on 
all three boundaries to new random values. The results 
for these boundary conditions give similar results to those 
for the R boundary condition. 

It is important to estimate the size of the errors which 
arise because the algorithm is not guaranteed to find the 
exact ground state. This is a problem only for the larger 
sizes so we estimated the errors carefully for i = 8 and 10 
for P-AP boundary conditions by doing a smaller num- 
ber of samples for a larger number of runs, n^ = 10, than 
before (see Table I) and assuming that the difference in 



results is a reasonable measure of the error in the orig- 
inal data. For L = 8, the result for whether the block 
spin configuration changed upon changing the boundary 
conditions was the same after 2 runs as after 10 runs in 
367 out of the 370 samples considered, with 1 sample giv- 
ing \q\ = 1 after 2 runs but \q\ < 1 after 10 runs, and 2 
samples the other way round. Hence our estimate of the 
relative error in the data in Fig. ^ (due to not finding the 
correct ground state) is in the range (— 1 ± \/5)/370, i.e. 
between —0.7% and -1-0.2%, which is well within the sta- 
tistical error bar of ±1.5% shown in Fig. O. For L = 10 
the corresponding figures are: 11 out of the 266 samples 
considered gave |g| = 1 after 1 run but |g| < 1 after 10 
runs, and 13 samples the other way round. The estimate 
of the error is therefore (— 2±\/24)/266, i.e. the true an- 
swer should lie between —2.5% and -1-1.1% of the value 
given, compared with the statistical error bar of ±2.5% in 
Fig. O. The error from the algorithm is therefore no big- 
ger than the statistical error, and is significantly less than 
that on the positive side. Hence the decreasing trend in 
the 3d data in Fig. y, is not due to inaccuracies in the 
algorithm. Altogether, we find no evidence for a system- 
atic error in the data in Fig. g due to not always finding 
the true ground state. A similar analysis for the defect 
energy with P-AP boundary conditions finds that, for 
L = 8, the error is negligible compared with the statis- 
tical error, but that for L = 10 the result is about 5% 
too high, compared with the error bar in Fig. |] of about 
±1.5%. This is probably why the L = 10 data point lies 
above the fit. The best fit with the L = 10 point 5% 
lower has 6 = 0.21, rather than 0.23, which is why we 
give asymmetric error bars in Eq. (|6|). 

In d = 2, the data for P-AP and P-R lie nicely parallel 
to each other and we find A = 0.69 ± 0.02, and hence 
df = 1.31 ±0.02, for both sets of boundary conditions, in 
agreement with PY who get A = 0.70 ± 0.08 and also in 
agreement with Middletor£3. Note that A is substantially 
smaller in d = 3 than in d = 2. It would be interesting to 
know how it varies in higher dimensions, but the numerics 
of such a calculation are challenging. 

The only other calculation oi df in d = 3 that we are 
aware of is that of HuseE3 who studied domain growth at 
finite temperatures by Monte Carlo simulations, obtain- 
ing df ~ 2.2. However, his method "only measures the 
size of domains which are compact" , which may explain 
why his value is lower than ours. 

To conclude we have seen that the ground state struc- 
ture appears to be trivial in a spin glass model with a 
finite Tc, the three-dimensional Ising spin glass. It re- 
mains to understand why Monte Carlo simulations at 
finite temperature find, by contrast, evidence for a non- 
trivial pure state structure. 
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